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The effect of non-universal SUSY soft breaking on predictions of dark matter de- 
tector event rates are surveyed for supergravity models with gravity mediated soft 
breaking. For universal soft breaking in the first two generations and tan fi ^ 20, 
non-universal effects can be characterized by four parameters (two for Higgs and 
two for third generation squarks) in addition to those of the minimal model 
(MSGM). These can increase or decrease event rates by a factor of 10-100 in the 
domain m^o ^ 60 GeV (x? = lightest neutralino) but produce generally small 
effects at higher masses. The value of the top mass and 6^5 + 7 branching 
ratio eliminates most of the parameter space for fj, < 0, causing event rates for 
\i < to be a factor 100 smaller than for fi > 0. A correlation between large 
(small) event rates and small (large) b — > s + 7 branching ratio is observed. The 
effect of future satellite (MAP and PLANCK) precision determinations of cosmo- 
logical parameters on predicted event rates is examined for examples of the ACDM 
and vCDM models. It is seen that these could sharpen the event rate predictions 
and also restrict the allowed gaugino mass ranges, thus influencing predictions for 
accelerator SUSY searches. 

1 Introduction 

The composition of the dark matter (DM) in the universe, which makes up 
90% or more of the universe's total matter, is one of the most important 
unresolved problems in astronomy. At present, dark matter has been observed 
only from its_gravitational interactions and may consist of a number of different 
components El. There may be baryonic (B) dark matter (e.g. machos), hot dark 
matter (HDM) which was relativistic at the time galaxy formation (possibly 
massive neutrinos) and cold dark matter (CDM) which was non-relativistic 
during galaxy formation. In addition a cosmological constant (A) may be 
present. 

The amount of each type of DM can be specified by the parameter £1^ = 
Pi/Pc where pi is the matter density of type "i", p c = 3H 2 /8itGn is the 
critical matter density to close the universe, H is the Hubble constant, H = 
hlOOkm/sMpc and Gn the Newtonian constant. Currently, measurements of 
h fall in the range 

0.5 < h < 0.75 (1) 
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and different models give estimates for the cosmological CDM of 

0.1 < ttcDAih 2 < 0.4 (2) 

The CDM that can be directly observed is that which exists locally in the Milky 
Way. This has been estimated to have a density of ~ 0.3 GeV/cm 3 , 

impinging on the Solar System with velocity Vdm — 300fcm/s. 

The cold dark matter is of particular interest in that if it is of particle 
nature it is most likely an "exotic" particle, i.e. one not found in the Standard 
Model (SM). We consider here models of physics beyond the Standard Model 
based on supergravity grand unification with R-parity invariance, where su- 
persymmetry (SUSY) breaking takes place at a scale near or above the GUT 
scale Mq — 1.5 x 10 16 GeV. This breaking occurs in a "hidden" sector and 
is transmitted (super) gravitationally to the physical sectorH. Such models au- 
tomatically predict the existance of CDM in that the lightest supersymmetric 
particle (LSP) is absolutely stable and over most of the parameter space is the 
lightest neutralino, Xv Thus the relic Xi left over from the Big Bang would 
be the CDM seen today. Further, over a significant part of the SUSY param- 
eter space, the amount of CDM predicted is in accord with what is observed 
astronomically, i.e. Eq.(2). 

2 DM Detector Event Rates 

Terrestial experiments detect incident local (Milky Way) DM particles by their 
scattering by quarks in nuclear targets. We briefly review in this section the 
relevant formulae for prediction of detector event rates for SUSY models. 

The analysis proceeds as follows. One first calculates the relic density 
of neutralino CDM n which remain after neutralino annihilation in the early 
universe. One finds Em: 

a «tf - 2.48 x 10-" (^V(— V^lI (3) 
\Tj \2.73j J(x f ) 1 ' 

where Xf = kTf/m x o (Tf = neutralino freezeout temperature, Nf = number 
of degrees of freedom at freezeout), T 7 is the cosmic microwave background 
(CMB) temperature, (T x o/T 7 ) 3 is the reheating factor and 

r x f 

J(x f ) = I dx(av)(x) GeV~ 2 (4) 
Jo 

In Eq.(4), a is the neutralino annihilation cross section (calculated from the 
SUSY model), v the relative velocity and ( ) means thermal average. One 
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sees that £l x oh 2 scales inversely with the annihilation cross section, i.e. the 
more annihilation there is, the less relic Xi remain. One restricts the SUSY 
parameter space so that Eq.(2) is obeyed, and also that the current SUSY 
parameter space bounds from LEP, the Tevatron and CLEO (6 — > s + 7 decay) 
are satisfied. 

One then calculates, in the restricted SUSY parameter space, the quan- 
tity R, the expected event rate of detector scattering events (pet kilogram of 
detector per day) for the incident flux of Milky Way neutralinosEl: 



R — (Rsi + Rsd) 



0.3 GeV/cm 3 



320/cto/s 



events 
kg d 



(5) 



where 



!6m x oM N M% 
Rsi = T77 — ~, 10" I Asi 



[M N + TO o] 



(6) 



Rsd = 



16to v oMn 
-*i 

[M N + m of 



r A 2 J(J 4- 1) I A S d | 2 



(7) 



where Mm is the nuclear target mass and J is its spin, Mz is the Z boson 
mass and (N\ £ S*i\N) = X(N\ J\N) with Si the ith nucleon's spin. In Eqs.(6) 
and (7), Asi(Asd) are the spin independent (spin dependent) scattering am- 
plitudes. Note that for heavy nuclear targets, Rsi ~ Mjv while Rsd ~ ^-/Mn 
making the heavier targets generally more sensitive than lighter ones and Rsi 
generally the dominant contribution for heavy targets. There are a number 
of uncertainties in the above analysis involving the strange quark contribution 
to the nucleon, the nature of nuclear form factors etc., making the theoretical 
predictions of R uncertain to perhaps a factor w 2. 



3 Soft Breaking: Universal Case 

We consider here models where the GUT group G breaks to the Standard 
Model group at M G : G -» SU(3) x SU(2) x [7(1). The simplest supergravity 
model of this type, the minimal SUGRA model (MSGM) depends at Mq on 
four extra soft breaking parameters and one..sign to determine all the masses 
and interactions of the 32 SUSY particles Era. These new parameters are mo 
(the universal scalar soft breaking mass), m 1 / 2 (the universal gaugino mass), 
Aq (the universal cubic soft breaking parameter), Bq (the quadratic soft break- 
ing parameter) and the sign of /zo, the Higgs mixing parameter in the effective 
potential term ^qH\H2 (where Hi,i = 1,2 are the two Higgs doublets). The 
renormalization group equations (RGE) then allow one to proceed downward 
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to lower energy scales where the soft breaking parameters trigger the break- 
ing of SU(2) x £7(1) El. In fact one may show that a necessary condition that 
electroweak breaking occur at a lower scale is that at least one soft breaking 
parameter and /Uo be non-zero at Mq, and this breaking will occur at the elec- 
troweak scale (Mz) provided m* obeys 90 GeV t nit t 200 GeV. Thus the 
model automatically requires a heavy top quark and it is SUSY soft breaking 
at Mq that gives rise to electroweak breaking at Mz- 
The conditions for electroweak breaking are 



2 - M2 tan 2 /3 1 2 . 2 2Bu, 

1 = 7^ ~M?: sin p = „ 7j — 

tan/3 2 -1 2 Z: ' \i\ + \i\ + 2/i 2 



where u 2 , /i 2 , fl are running parameters at Q = M^, /j, 2 = m 2 H . tan /3 = 

(H2) /(Hi) and ^ are 1°°P corrections Thus radiative breaking allows one 
to determine /it 2 and eliminate B in terms of tan /3. The determination of | \x \ 
greatly enhances the predictive power of the model. In the following we will 
restrict the parameter space examined to be in the domain 

m , mg < ITeV, 2 < tan/3 < 25, -7 < A t /m < 7 (9) 

where rrig is the gluino mass and A t is the i-quark A-parameter at the elec- 
troweak scale. (Eq.(9) satisfies usual "naturalness" conditions.) For_most of 
the parameter space, /i 2 /M| 3> 1 which leads to "scaling" relations for the 
charginos (xf^ = 1, 2) and neutralinos z = 1, .., 4): 



2m x o = m x ± = rn x o = ( - - -)m g ; m g = m 1/2 (10) 



m ± = m^o = m^o 3> m^o (11) 

X 2 a 3 X 4 Xi v / 

In addition, the four Higgs bosons, h, H°, A° and i? ± obey in this domain the 
ralations 

m H o = m A = m H ± » m h (12) 

and m h < 120 GeU. 

Most of the SUGRA dark matteLanalysis has been done within the above 
framework of universal soft breakingcra. We examine next what modifications 
arise when non-universal soft breaking occurs. 



4 Soft Breaking: Non-Universal Case 

The MSGM model with universal soft breaking parameters, is of course the 
simplest SUGRA model, which in part is why it has been examined so exten- 
sively. However, there are a number of reasons why one might expect non- 
universal SUSY soft breaking to occur at Mg- Thus in general, non-universal 
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soft breaking will arise if in the Kahler potential, the interactions between the 
hidden sector fields (whose VEYs give rise to SUSY breaking) and the physical 
sector fields are not universal Further, even if universality where to hold at 
a more fundamental level, e.g. at the Planck or string scales, running the RGE 
dowzLto the GUT scale can produce significant non-universal soft breaking at 
MgEJ. Finally, we note that in the breaking of a higher rank GUT group down 
to the SM group at Mq, the D terms can generate non- universalities ElL 

Flavor changing neutral currents [FCNC] will be suppressed iLthe soft 
breaking masses of the first two generations are universal at Mq Ej, and in 
the following we will assume that this is the case and let mo be their common 
mass. We then parametrize the Higgs and third generation squark and slepton 
masses at Mq as follows: 



(13) 



m q L = m o(! +^3); 



m 2 (l + <5 4 ); < =m 2 (l + S 5 ) (14) 



(15) 



where — {ul,<Il) is the left squark doublet, li — (vl^l) the left slepton 
doublet etc. The Si represent the deviations from universality, and we will 
restrict these to obey — 1 < Si < +1. In addition, we denote the A parameters 
at Mq by Ao t , Aqj, and Aq t which also need not be universal. 

For GUT groups containing an SU (5) subgroup with matter embedded in 
the 10 + 5 of SU{5) (e.g. SU(N), N > 5; SO(N), N > 10; E 6 ) one has 



63 = 64 = 65; 6 6 = S 7 ; A ob = A 0t 



(16) 



We note that for tan/3 ^ 20 , 6 5 ,6 6 ,6 7 and A b,A 0T do not enter significantly 
in the calculations, and we will set them to zero in the following. 

To examine the significance of the non-universal soft itreaking, we exhibit 
the 1-loop RGE expressions for /x 2 of Eq.(8) for Q = M z 



t 2 



1-3D, 
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(17) 
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where t = tan/3, D = 1 - (m t /200 sin/3) 2 , A R = A t - 0.61m 5 , a G ^ 1/24, 
5o = TrYm 2 (Y = hypercharge and m are the soft breaking masses at Mg) 
and = §D (1 - Do){HJF) 2 + e- g/t 2 (with the RGE form factors H, F, 
e, g given in Ibahez et aid). In Eq.(17), the m 2 , term has been divided into 
a universal part and a non- universal part (dependent on the Si). D vanishes 
at the i-quark Landau pole and is generally small (0 < Do ^ 0.25). An is 
the residue at the Landau pole (i.e. A t = Ar/Do — (H 3 /F)m 1 / 2 )- Since 
nii/2 = (ac / Ofijnig, where ac is the GUT coupling constant, the C^m 2 ^ 
term scales with m- or alternately by Eq.(10) with (m x o) 2 . So vanishes for 
universal soft breaking masses (by hypercharge anomaly cancelation) but is 
non-zero for non-universal masses. (The So term is generally small, usually 
only a few percent correction.) 

5 Effects Of Non-Universal Soft Breaking On DM Event Rates 

While many parameters enter into the DM event rate formula Eq.(5), /i 2 plays 
a central role in determining R. This is because /i 2 governs the interference 
between the Higgsino and gaugino parts of \i m the x? " quark scattering 
amplitude and it is this interference which gives rise to Rsi (which we have 
seen for most detectors is the dominant part of R). One finds in general that 
increasing (decreasing) decreases (increases) the amount of interference and 
hence the size of R. 

Since Do is small, one sees from Eq.(17) that S 3 and £4 contribute oppo- 
sitely to 82 in 11 2 . Thus in evaluating the contribution of non-universal masses, 
it is necessary to consider both Higgs and squark masses since they produce 
effects of the same size. There are certain situations where the non-universal 
terms become relatively large compared to the universal contributions, greatly 
enhancing the non-universal effects. Thus the universal contribution of the 
Landau pole term (which is generally quite large) vanishes if the residue An 
is zero, i.e. when A t = 0.61wg. Also the term is small for light \i 

(or light g) according to the scaling relations Eq.(10). In addition, for t 2 ^ 1 
(i.e. t ^ 3), the universal contribution to the m 2 , term becomes small when 
D w 1/3. 

The above ideas allow one tn. understand qualitatively the detailed com- 
puter calculations of event ratesEij. Fig.l shows maximum and minimum event 
rates for the case 63 = = 64. The solid curves are for universal soft break- 
ing. From Eq.(17), the case 82 = — 1 = —8% (dotted curves) corresponds to 
increasing /i 2 , which reduces the event rates. This reduction can be as much 
as a factor of O(10 — 100) for m x o ^ 60 GeV . The effect is largest for the 
minimum event rates, since these occur at small tan (3 (which by Eq.(17) would 
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40 60 80 100 120 

Neutralino Mass in GeV 

Figure 1: Maximum and minimum event rates for a Xe detector as a function of fn a for 
At > 0, 0.1 < n x oh 2 < 0.4 with <5 3 = = <5 4 and <5 2 = = <5i (solid), <5 2 = -1 = 1 -(5 i 
(dotted), <5 2 = 1 = -<5i (dashed) 0. 

magnify the effect). The reverse situation, 82 = 1 = — Si (dashed curve) causes 
a decrease of /z 2 and hence increase of event rates. The effects of non-universal 
soft breaking masses becomes small for heavy neutralinos, m x o > 60 GeV, 
since then the C fl m 1 / 2 becomes larger masking the non- universal effects. 

Fig. 2 shows a similar analysis with Si, 82 and 83 set to zero. We see the 
curves with 64 = 1 (dotted) resembles the 62 = — 1 = —S\ curves of Fig. f while 
the 84 = —1 (dashed) resembles the 82 = 1 = — #1 curves of Fig.l. Fig. 3 shows 
curves similar to Fig. 2 with 83 non-zero. These results follow from Eq.(17) 
where the effects of S3 and <5 4 were seen to be opposite to those of 8 2 , e.g. a 
negative £4 simulating a positive 82 etc. 

As seen in Eq.(f6), 83 = 84 for most GUT groups. Fig. 4 examines this 
case with the choice £3 = 84 = — 1 . Here the squark non-universal masses add 
coherently in Eq.(f7) with the Higgs non- universal masses for the case 82 = 
1 = —81 to reduce /x 2 . The maximum event rate for this case (dashed curve) 
which occurs for large tan/3, is then greatly enhanced over a wide range of m x o 
rising to R ~ £_L— 10) events / kg d. This value is the current sensitivity of the 
Nal detectors EJ, and so this case is at the edge of experimental observation. 
For the reverse situation, 82 = — 1 = —Si (dotted curve), the squark and Higgs 
non-universal effects mostly cancel in Eq.(17) and one indeed finds that these 
curves lie close to the universal one for most of the neutralino mass range, even 
though a large amount of non-universal soft breaking has occured. 
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40 60 80 100 120 

Neutralino Mass in GeV 



Figure 2^_Same as Fig.l with <5i = 82 = £3 = and 64, = (solid), 64 = +1 (dotted), S = — 1 
(dashed) bl. 




40 60 80 100 120 

Fig. 3. Neutralino Mass in GeV, O.KOmega h-2<0.4 



Figure 3: Same as Fig.l with <5i 
«5 3 = -1 (dashed). 



and <5 3 = (solid), <5 3 = +1 (dotted), 
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Neutralino Mass in GeV 

Figure 4: Same as Fig.l with 8x_= — 1 = <?4 and S2 = = 8± (solid), 82 = —I = —81 
(dotted), 8 2 = 1 = -Si (dashed )B 



The above discussion shows how squark and Higgs non-universal soft break- 
ing masses can interact with each other, sometimes enhancing and sometimes 
canceling their effects depending on the sign of the non-universal terms. 

6 The t Quark Mass, b -> s + 7 Decay And The Sign Of fx 

While the value of /x 2 is a crucial parameter in determining DM event rates, 
non-universal soft breaking masses also enter in other parameters which control 
the event rates in a more indirect way. Thus the stop mass matrix is given by 

M 2 _ ( m 2 tL -mt(A t + fJ, cot P)\ , , 

M t ~ {- mt (A t + fi cot 0) m \ R )■ 

where 

m t L = m Q L + m t + ( 5 _ 3 sin2 Q w)M 2 z cos 2/3 (19) 

2 

m 2 tR =mfj+m 2 + - sin 2 ^> W M\ cos 2/3 (20) 
o 

Here El 

-kl-n )§ + C Q rni /2 -hl-^l)S (21) 
9 



- D Q )^ + Cam? /2 + 1(1 - ^M)5 (22) 

and Cq, Cjj are given in Refill. Note here S3 and 64 act oppositely to each 
other, though if the GUT condition 83 = 84 is imposed, they still act oppositely 
to S 2 . 

Since the top quark is heavy, i.e. Do is small, the Landau pole term in 
Eq.(22) can drive fi, theJight stop, tachyonic if Ar is large enough, eliminating 
such parameter points uJ. Since A 2 R = (A t — O.dlrrig) 2 , regions where A t is 
negative get eliminated, and one finds the allowed region of parameter space 
is restricted byEBIij 

A t /mo t -0.5 (23) 

Thfipii — > s + 7 decay also strongly restricts the SUSY parameter-, 
space Oo The measured branching ratio for this process by CLEO is 113 
B(B -> X sl ) = (2.32 ± 0.67) x 10~ 4 . The SM prediction, with NLO effects 
included is B(B — > X s j) = (3.48 ±0.31) x 10 -4 where the error in the theoret- 
ical calculations is now dominated by uncertainties in the input parameters, 
e.g. a$(Mz), nit, m c /mb, etcEa While these two results are statistically con- 
sistent, it is clear that there is a certain amount of tension between the data 
and the SM. The SUSY corrections to the SM reduce the SM branching ratio 
for A t / n > and increase it or A t / [i < 0. Since by Eq.(23) A t is mostly posi- 
tive, and the experimental branching ratio lies below the SM value, the current 
CLEO data already eliminates most of the /i < part of the parameter space 
at the 95% C.L. level. The part of the u, < parameter space that survives is 
for small tan /?, and so one finds in most of the remaining parameter space, the 
event rates for /i < are a factor of 100 (or more) smaller than those of /i^ 
and this result is true with or without non- universal soft breaking masses E3. 

A second feature involving the b — > s + 7 decay is the existance of a 
correlation between the DM event rate R and the b — * s + 7 branching ratio 
B: large (small) event rates tend to correlate with small (large) branching 
ratio. This shows up most strongly if one considers the maximim (R max ) or 
minimum (R m i n ) event rate at fixed values of the neutralino mass. One finds 
(for universal soft breaking) that the part of the parameter space with B < 
2 xlO -4 corresponds to R ma x ~ 0.1 events/kg d, and B > 3x 10~ 4 corresponds 
to R m ax ~ 0.01 events/kg d. Similarly for R m in < 0.003 events/kg c?, one 
finds B > 3 x 10~ 4 . The correlation-is less strong for values of R in between 
Rmax and Rmin, but still significant Eil. 
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It is clear from the above that as the b — > s + 7 data becomes more 
precise, it will have a strong impact on DM predictions, and on restricting the 
supersymmetry parameter space in general. 



7 Determination Of Cosmological Parameters 

One of the major sources of uncertainty in dark matter event rate predictions 
is the lack of accurateJinowledge of the basic cosmological parameters. Future 
satellite experiments E3, MAP (scheduled for the year 2000) and PLANCK 
(planned for 2005) will be able to measure the angular power spectrum very 
accurately. Thus one defines the correlation function 

Jo Jo 

where Tq — 2.728 ± 0.002°AT is the cosmic microwave background (CMB) 
temperature, AT(qi), i = 1,2 are the deviations from To in directions qi and 
cos 9 = q± ■ <72 ■ Expanding in Legendre polynomials 

0/ 1 1 

c(#) = V— — qP/(cos0) (25) 

* — ' 47T 

these experiments will be able to determine the q out to I w 1000. Since the 
q are sensitive to the different cosmological parameters, it will be possible to 



determine the Hubble constant, the amount,, 
constant etc. to an accuracy of a few percentt 



ark matter, the cosmological 
. We consider here what such 



determinations might mean for DM predictions for two cosmological models. 



7.1 A CDM Model 

Current astronomical measurements suggest that while there is a large amount 
of CDM, it's mean density is considerably less than the critical density p c 
and that models with a cosmological constant A and Sic dm ~ 0.4 are good, 
if not better fits to the current astronomical data E£l. We consider here the 
ACDM model with cold and baryonic DM such that Sltotal = 1 where tltotai = 
dm + + f^A- We assume that the central values of the cosmological 
parameters measured by PLANCK and MAP are 

Q CDM =0A, n B =0.05, fi A = 0.55, /i = 0.62 (26) 

(consistent with current determinations of these quantities) . The accuracy with 
which each of these quantities can be measured by the-PLANCK satellites has 
been estimated in RcfE3 and one finds from this that til 

n CDM h 2 = 0.154 ±0.017 (27) 
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ACDM, 5 =0=5, 
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300 400 500 600 

Gluino mass (GeV) ; n >0 

Figure 5: Maximum and minimum event rates for a Xe detector for /i > as a function 
of rrig for universal soft breaking. The solid (dashed) curves are the lstd (2std) range of 
Eq.(27)0. 

This is to be compared with the very broad window of Eq. (2) for the current 
estimates of Hcdm^ 2 - 

Assuming as before that the CDM is the relic neutralinos, this sharpening 
of the allowed region of tt x o h 2 significantly effects the DM event rate predic- 
tions. Fig. 5 shows the maximum and minimum event rates as a function of rrig 
for universal soft breaking (Si — 0). (Recall m x a scales with rrig by Eq.(10).) 
Comparing with e.g. the solid curve of Fig.l, we see that the effect of Eq.(27) 
is to increase R m in significantly and put an upper bound on rrig (and hence 
on m x o ). This latter effect arises from the fact that the \i annihilation cross 
section in the early universe is a decreasing function of m x o, and hence an 
upper bound on the relic density produces an upper bound on m x o. Fig. 6 
shows the effects of non-universal soft breaking on the event rates. Note that 
for 62 = 1 = — Si there is also a minimum value of m g (and hence a minimum 
rn v o ) since for this case a lower value of m^o would violate the lower bound on 
Q x oh 2 of Eq.(27). (The additional isolated points for 5 2 = 1 = —Si arise from 
the fact that for this case /i 2 is driven small (as can be seen from Eq.(17)) and 
so the scaling relations Eq.(10) no longer hold allowing m x o to be sufficiently 
small so that the upper bound of Eq.(27) can be satisfied even though rrig is 
large.) We also note that for S2 = 1 = — Si, almost the entire parameter space 
has now been eliminated for /i < 0. 
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Figure 6: Same as Fig. 5 for 1 std range of Eq.(27) with 82 = = <5i (solid), 82 = —1 = — <5i 
(dotted), 82 = 1 = —Si (dashed) and 83 = = 84. The isolated points are for 82 = 1 = — Si 
(Rmax are empty squares, R m i n are solid squares)o 



The above discussion shows that the upper bound on fi x o/i 2 of Eq.(27) for 
the ACDM model produces upper bounds on m g and by scaling, on the other 
gauginos. Thus (aside from the discrete points in Fig. 6) the lstd (2std) gaug- 
ino upper bounds are m g < 520(560) GeV, m x o < 70(77) GeV, m x ± t 150 

GeV, and also the light Higgs obeys rrih t 120 GeV. It is interesting to com- 
pare these bounds with what might be the expected reach of an upgraded 
Tevatron with 25/6~ 1 of data. Thus it is estimated that the gluino could 
be observed with a mass reach of rrig * 450 GeV for most of the parameter 
space, the chargino with m^± * 150 GeV for about 2/3 of the parameter space 

and a Higgs for rrih t 120 Gey@. Hence the example of the ACDM model 
considered here would suggest that an upgraded Tevatron could examine a sig- 
nificant part of the SUSY parameter space allowed by the PLANCK satellite 
measurements of the cosmological parameters. 

As discussed in Sec. 6, there is a correlation between large (small) DM event 
rates R and small (large) b — > s + 7 branching ratios B. This is exhibited 
in detail for the ACDM model in Fig. 7 which shows a general scatter plot 
over the full parameter space of the model of Eq.(27). Recall that the curent 
CLEO data is B = (2.32 ± 0.67) x 10" 4 while the SM prediction is EJ 
B = (3.48 ±0.31) x 10~ 4 . One sees from Fig.7 that if B > 3 x 10~ 4 (consistent 
with the SM predictions) then R is small i.e. R < 0.1 events/ kg d while if 
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Figure 7: Scatter plot over allowed SUSY parameter space for R vs. B (b — » s + 7) for 
ACDM model for fx > with universal soft breaking masses and lstd bounds of Eq.(27). Ed 

B < 3x (suggested by the current CLEO data) then R > 0.05 events/kg d 
for almost all parameter points. Thus one finds an interesting correlation 
between accelerator physics and cosmology. 

7.2 vCDM model 

If neutrinos do have masses in the eV range, they would contribute significantly 
to the hot dark matter of the universe. As a second example of what the new 
satellite experiments might see, we consider the neutrino cold dark matter 
(j/CDM) model and assume that the MAP and PLANCK satellites determine 
the following central values for the cosmological parameters: 

0, = 0.2; ftcDM = 0.75; n B = 0.05; h = 0.62 (28) 

Using the estimateicLaccuracies that these quantities can be measured by the 
PLANCK satellite BEfJ one finds Ed 

QcDMh 2 = 0.288 ±0.013 (29) 

As in the ACDM model, one finds a narrowing in the difference between the 

maximum and minimum event rates for m v o 'Z 60 GeV . There is also an 

1 

upper bound on the gaugino masses, which now is higher than in the ACDM 
model since VL x oh 2 is larger. One finds for the lstd (2std) bounds of Eq.(29) 
thatE3mg t 700(720) GeV, m x o < 95(100) GeV and m x ± < 200 GeV. Thus 



14 



for this model, the LHC would be crucial in order to explore the full SUSY 
parameter space. The lstd window of Eq.(29) also shows an additional feature 
of gaps in the allowed range of m g (and m x o). While these gaps fill in at 
the 2std level, they illustrate the impact that astronomical measurements may 
have on accelerator physics searches. 

8 Conclusions 

We have considered here the effects of non-universal SUSY soft breaking on 
predictions of dark matter detection rates for supergravity models with R 
parity where supersymmetry is broken at a scale Z Mq, the breaking being 
communicated by gravity to the physical sector. The lightest neutralino be- 
comes the CDM candidate for almost the entire parameter space. For models 
which have universal soft breaking in the first two generations (to suppress 
FCNC processes), one generally needs nine additional parameters to describe 
the non-universal effects of the Higgs boson and third generation squark and 
slepton masses (which reduce to five parameters if one imposes the symmetry 
constraints, Eq.(16), of most GUT groups). For tan/3 < 20 only four of these, 
two Higgs masses and two squark masses, enter significantly. While this is 
still a major enlargement of the minimal SUSY parameter space, the extensive 
knowledge already obtained from studying SUGRA dark matter predictions 
for universal soft breaking have allowed one to determine in large measure the 
effects that such non-universal soft breaking produce. Thus the squark and 
Higgs non-universal effects can act to enhance or cancel each other (depending 
on their relative signs) and can increase or decrease event rates by a factor of 
ps 10 - 100 in the domain m v o < 60 GeV (ras < 400 GeV ). However, the 
non-universal effects are generally small at higher masses t3. 

SUGRA models can account for both the existance of CDM as well as 
accelerator SUSY phenomena, and interaction between the two sets of phe- 
nomena has begun to occur. Thus the fact that the top quark is heavy, and 
that the measured b — ► s + 7 branching ratio lies about 1.6 std below the SM 
prediction eliminates most of the \x < and A t < parts of the parameter 
space at the 95% C.L. 15-18 This then leads to DM event rate predictions for 
the small remaining part of the parameter space with \x < to be ps 100 times 
smaller than for /i > 0E3. Further one finds a correlation between large (small) 
event rates and small (large) predicted b — ► s + 7 branching ratio for SUGRA 
modelscd. One expects that further b — > s + 7 data will play an important role 
in SUGRA dark matter event rate predictions. 

Future satellite experiments by PLANCK and MAP as well as ballon and 
ground based experiments will be able to determine the basic cosmological 
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parameters with great accuracy (at the 1-10% level) and hence greatly reduce 
the uncertainties in DM predictions. To illustrate what might be expected 
from the PLANCK and MAP experiments, examples of the ACDM model 
and z^CDM model were considered. The narrowing of the allowed window 
for fl x oh 2 generally reduces the spread between the maximum and minimum 
event rates for m v o ^ 60 GeV. and limits the maximum (and for some signs of 
non-universal soft breaking the minimum) values of gaugino masses EiJ. Thus 
astronomical measurements can have a significant impact on SUSY accelerator 
searches. 
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